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Abstract 

We numerically study the effects of varying electric conductivity and magnetic permeability of 
the bounding wall on a kinematic dynamo in a sphere for parameters relevant to Madison plasma 
dynamo experiment (MPDX). The dynamo is excited by a laminar, axisymmetric flow of von 
Karman type. The flow is obtained as a solution to the Navier-Stokes equation for an isothermal 
fluid with a velocity profile specified at the sphere's boundary. The properties of the wall are taken 
into account as thin-wall boundary conditions imposed on the magnetic field. It is found that an 
increase in the permeability of the wall reduces the critical magnetic Reynolds number Rm cr . An 
increase in the conductivity of the wall leaves Rm cr unaffected, but reduces the dynamo growth 
rate. 
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Over the past decade, significant effort has been directed at the experimental demonstra- 
tion of dynamo action - self-excitation and maintenance of the magnetic field in a flowing 
electrically conducting fluid. A number of experiments with liquid metals have been con- 
structed to test this phenomenon in various settings [1-6], and successful observations of 
dynamo action have been reported in three of them [4-6]. Among other things, these ex- 
periments revealed the critical importance of the magnetic properties of the flow-driving 
impellers. Namely, the von Karman sodium experiment only self-sustained a dynamo field if 
the impellers were ferromagnetic [6, 7]. In addition, the finite resistivity of the experimental 
container is expected to be crucial for the dynamo instability - a situation similar to the re- 
sistive wall mode (RWM) in tokamaks [8]. Normally stable for the perfectly conducting wall, 
the RWM can become unstable if the wall has finite resistivity; in this case the instability 
develops on the wall's resistive time scale. These facts initiated more thorough theoretical 
studies of the effect of the imposed boundary conditions on the dynamo in experimentally 
relevant models [9-16]. The studies established that there are no general dependences of 
dynamo properties on conductivity and permeability of the boundary; the dependences are 
different for different models and flows. 

This circumstance motivates us to perform an analogous study for the Madison plasma 
dynamo experiment (MPDX, Fig. 1), currently under construction at the University of 
Wisconsin-Madison. The experiment is aimed at investigations of fundamental properties 
of dynamos excited by controllable flows of plasmas. Its original design was proposed in 
Refs. [17, 18] and conceptual features were successfully tested in the plasma Couette ex- 
periment (PCX) [19]. The experimental vessel is a sphere of 3 meter in diameter. An 
axisymmetric multicusp magnetic field (created by 36 equally spaced rings of permanent 
magnets with alternating polarity) confines the plasma. The field is localized near the vessel 
wall and a large volume of unmagnetized plasma occupies the experiment's core. An electric 
field applied across the multicusp field drives the edge of the plasma azimuthally. Arbitrary 
profiles of azimuthal flow v^(0) can be imposed at the spherical boundary by modulating 
the electric field as a function of polar angle 9 using discrete electrodes. 

Results of Ref. [18] show that some flows generated in such a way can lead to a dynamo 
instability. However, in Ref. [18] insulating boundaries are assumed, whereas in MPDX the 
vessel is made of aluminum whose conductivity is much higher than that of the plasma 
under expected experimental conditions. The goal of this paper is to generalize the results 
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FIG. 1: Madison plasma dynamo experiment (MPDX): (a) sketch of the experiment; (b) electrode 
configuration near the wall for driving plasma velocity v^(0); (c) model for deriving thin- wall 
boundary conditions, with a and fi denoting conductivity and relative permeability of the respective 
media. 

of Ref. [18] by considering effects of varying conductivity and permeability of the vessel on 
the dynamo in a model relevant to MPDX. 

To describe the plasma we use dimensionless numbers: 

M = V M, Re=^ Rm=^,Pm=- 
V Po v V V 

- Mach, fluid Reynolds, magnetic Reynolds and magnetic Prandtl, respectively. Here Vo 
is the peak driving velocity, p and P are the average plasma mass density and pressure, 
Ro is the radius of the sphere (a unit of length throughout the paper), v and 77 are the 
plasma kinematic viscosity and magnetic diffusivity (assumed to be constant and uniform). 
For given plasma parameters these numbers can be estimated from the Braginskii equations 
[20] (see corresponding formulas in Refs. [18, 21]). Their expected values for MPDX are 
listed in Table I. By varying temperature, density and ion species of the plasma one can 
change its magnetic Prandtl number by several orders of magnitude. Such flexibility makes 
it possible to demonstrate a dynamo in a laminar flow by choosing a regime with Pm ~ 1 
and Rm ~ Re ~ 10 2 . This is an advantage over the liquid metal dynamo experiments, 
where Pm ~ 10 -5 and the flows are always turbulent. 

Our first step is to find an equilibrium velocity field capable of dynamo action. For 
simplicity, we do not focus on the details of plasma driving near the wall. We neglect 
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TABLE I: Expected parameters of MPDX 



Quantity 


Symbol 


Value 


Unit 


Radius of sphere 




1.5 


m 


Wall thickness 


d 


0.05 


m 


Peak driving velocity 


V 


0-20 


km/s 


Average number density 


n 


10 i7 _ 1Q 19 


m 


Electron temperature 


T 


2-10 


eV 


Ion temperature 


T 


0.5-4 


eV 


Ion species 




H, He, Ne, Ar 




Ion mass 


Hi 


1, 4, 20, 40 


amu 


Mach 


M 


0-8 




Fluid Reynolds 


Re 


0- 10 5 




Magnetic Reynolds 


Rm 


- 2 x 10 3 




Magnetic Prandtl 


Pm 


10~ 3 - 5 x 10 3 





the multicusp magnetic field and applied electric field in our consideration and assume 
that the velocity profile is specified at the boundary. As shown in Ref. [21] for the model 
relevant to PCX (cylindrical prototype of MPDX), the velocity structure obtained under 
such assumption is the same as the velocity structure obtained with a more realistic E x B 
forcing, except in a thin boundary layer. 

The velocity field is found using the hydrodynamic part of the extended MHD code 
NIMROD [22] with an isothermal fluid model, which in non-dimensional form is 

HTI 

£T = -V-(nv), (1) 

„_ = _ n(v . V)v _ _ + _ (vV + jV(V • v) j , (2) 

where r, n and v stand for normalized time, density and velocity, respectively: r = t Vo/Ro, 
n = p/po, v = V/Vo- The differential plasma driving near the wall of MPDX is represented 
by the velocity boundary condition: 

v| r=1 = v^{6)e^ < 9 < 7T, (3) 

where v^iO) is a function of polar angle 6 with physical restriction ^(0) — v^(tt) — 0. In 
general, this function may be expressed as v^(0) = ^ a ^ sm k6 ) where a& are real coefficients. 
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(a) Velocity boundary condition 
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FIG. 2: Axisymmetric equilibrium flow of von Karman type for Mach number M — 1 and fluid 
Reynolds number Re — 300 used in kinematic dynamo study: (a) velocity boundary condition 
v $(0) adopted from Ref. [18]; (b) structure of normalized velocity; (c) contour plot of normalized 
density (dashed lines denote n < 1). Left half of (b) shows stream lines of poloidal flux nv po i 
superimposed on its absolute values depicted in colors, right half of (b) shows contour plot of 
azimuthal velocity (dashed lines denote < 0). Vertical lines in (b) and (c) represent the axis 
of symmetry. 

We use a velocity boundary condition of the von Karman type from Ref. [18], shown to 
result in a dynamo in an incompressible flow with Re = 300 and Rm > 237. It is given by 
a 2 = —0.4853, = —0.5235, a 6 = —0.0467, a 8 = 0.1516 (Fig. 2a). In present study we take 
the Mach number M = 1, the fluid Reynolds number Re = 300 and the magnetic Reynolds 
numbers up to Rm = 400. These parameters can be achieved in MPDX by creating, for 
example, an argon plasma with Vq = 5 km/s, no = 10 18 m -3 , T e = 10 eV and X L — 1 eV. 

In the NIMROD simulation, we used a meshing of the poloidal plane with 4608 quadri- 
lateral finite elements of polynomial degree 2, and 6 Fourier harmonics in the (^-direction 
(the azimuthal mode numbers are < m < 5). This resolution was sufficient for the laminar 
flow under consideration. We took a non-moving fluid (v = 0) with uniform density (n — 1) 
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as the initial state and evolved Eqs. (1), (2) with the boundary condition given by Eq. (3) 
until a steady state was reached. The resulting flow for Re = 300 and M — 1 is shown in 
Fig. 2. The velocity field v(r, 6) is axisymmetric and hydro dynamically stable with respect 
to perturbations with m > 0. 

The main results of the paper are obtained by solving the kinematic dynamo problem 
with this velocity field, 



7B = RmV x (v x B) + V 2 B, V • B = 0, 



(4) 



for unknown magnetic field B and normalized dynamo growth rate 7 = YR^/r]. We represent 
the divergence- free field as an expansion in a spherical harmonic basis [23] : 

1(1 + l)SiYr 



B ft 



B<h = 



E 

l=m 
L 

E 



E 

l=m 



Jmcf) 



1 dSi dY™ + imT x Y{ 



r dr 36 



r sin# 



imY™ dSi Ti dY r 



rsinO dr r 89 



D im(/) 



im<t> 



(5) 
(6) 
(7) 



where Si(r) and Ti(r) are functions of r only and Y^iO) are spherical harmonics related 
to the associated Legendre polynomials by Y l m (9) = P™(cos#). Since the velocity is ax- 
isymmetric, we consider each azimuthal mode m separately. The summation in Eqs. (5)- (7) 
is truncated at some spherical harmonic L (L = 20 provides a satisfactory convergence in 
these studies). Substituting Eqs. (5)- (7) into Eq. (4) and using the orthogonal properties of 
spherical harmonics, one obtains for m < I < L: 



r(2) 



8Sj 



+ 



r(3) T 

d 2 T x 1(1 + l)T t 
d 



J =171 



L 

-RmA^Y^ 

j=m 



(8) 



7 (i) T . 



dr 



r(2) T 



7-(3) 9Sj r(4) q \ , r(4) dSj 



(9) 



Here the bar above a symbol denotes its complex conjugate, Af 1 is a numerical factor 

(2/ + l)(/-m)! 



AT = 



2l(l + \)(l + m)V 



and iy 4 \r) are functions of r given by the integrals: 



r (2) 



r(3) 



r(4) 



j(i + 1) 
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U1 1 n ■ -\rr 

vq — sin 6/ — imv ( j ) Y l 



de 



d9, 



7T 

/ 



= «m 



7T 

J de Yl Y > M ' 



sin# 



de, 



j(j + 1) 



i 



v 4>~qq~ Sin ^ + imv e Y l 



de. 



(10) 

(11) 

(12) 
(13) 



Note that Eqs. (8)- (13) are valid for any axisymmetric velocity field. To calculate the 
integrals in Eqs. (10)-(13), we interpolate the velocity field on a uniform polar grid (typically 
with N r = 50 radial and Nq = 1000 angle grid points) and use the trapezoidal rule of 
integration. 

Eqs. (8), (9) should be supplemented with boundary conditions for functions Si(r) and 
Ti(r). The absence of a singularity in the field at the center of the sphere requires 



Si\ n = 0, TA =0. 

' I r=0 ' L \r=0 



(14) 



The outer boundary conditions depend on the properties of the shell. To avoid undesired 
diversion of flow-driving current into the shell, the inner surface in MPDX is covered with 
an insulating coating (Fig. lc). Thus, the normal component of current is zero at r = 1, i.e., 



\r=l 



0. 



(15) 



To derive the condition for Si at r — 1, we consider the model shown in Fig. lc and use the 
general boundary conditions for normal and tangential components of the magnetic field at 
the interface between two media with different relative magnetic permeabilities \i\ and /x 2 : 

D _ D B ^ _ B 2t 

& In — &2n) — • 

We also assume that the insulating coating is thin enough that it has no impact on profile 



7 



of Si. Then the resulting equations are (omitting in Si) 

ids i ds w 

r = l:S = S w , --£- = ^— , 16 

1 < r < 1 + — : —7^ = , 17 

Rq r] w or 1 r z 

r = l + — : S w = b v , — = , (18) 

Rq /jl w or or 

r > 1 + -J- : 5 V oc r" 2 , (19) 
Ro 

where Eq. (17) is derived for a stationary wall with thickness d, symbols with subscripts 
refer to wall (V) and vacuum ( 4 V), and symbols without subscript refer to plasma. We 
assume that d is small, so that the variations of Si in the wall are small too. This is the thin- 
wall approximation [14], it applies if d <C Rq and d <C l^/rl 1 / 2 . Under these assumptions, 
Eqs. (16)-(19) are reduced to 

^{l + lc^ + S{lfi + 7 c a )j =0, (20) 

where we have used the relation r\ = c 2 /(Air a /a) between magnetic diffusivity rj and electric 
conductivity a of a medium (c is the speed of light), and introduced the wall conductivity 
parameter c a and the wall permeability parameter c M , 



C ' = ^V c - = iv (21) 



Eq. (20) is obtained for a stationary wall without requiring the no-slip boundary condition 
for plasma velocity, in contrast to analogous equation (14a) from Ref. [14]. 

Eqs. (8), (9), (14), (15), (20) constitute an eigenvalue problem for the dynamo growth rate 
7 and unknown eigen- functions Si and T\. In order to solve it, we apply the finite difference 
method to Eqs. (8), (9) and discretize Si and 7} for each harmonic I (m < I < L) on a 
uniform grid at < r < 1 with N r equal intervals. Eqs. (14) and (15) are straightforward to 
implement in the finite difference scheme, while Eq. (20) is taken into account by using an 
extra (ghost) grid point to approximate the derivative. The resulting system of (L — m + 
l)(2N r — 1) linear algebraic equations for (L — m + l)(27V r — 1) unknowns is cast in the 
form of a matrix eigenvalue equation, which is solved in MATLAB. The developed scheme 
has been successfully benchmarked against the results of the kinematic dynamo study from 
Ref. [24]. 
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FIG. 3: Dependence of dynamo growth rate 7 on magnetic Reynolds number Rm for different 
values of the wall parameters c a and c M . 
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FIG. 4: Dependence of critical magnetic Reynolds number Rm cr on the wall permeability parameter 

Here we report the results of solving the kinematic dynamo eigenvalue problem [Eqs. (8), 
(9), (14), (15), (20)] with the velocity shown in Fig. 2, and for the relative permeability of 
the plasma \i — 1, magnetic Reynolds numbers Rm = — 400 and varying wall parameters 
c a and c M . We consider only the most unstable (or least decaying) m = 1 azimuthal mode. 
The convergence of the numerical scheme is checked by comparing simulations at different 
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resolutions. Results reported here are obtained by using a maximum number of spherical 
harmonics L = 20 and number of radial grid points = 50. 

The results are summarized in Figs. 3 and 4. In the present case 7 is always real, so 
the dynamo threshold Rm cr corresponds to the condition 7 = 0. As seen in Fig. 3, Rm cr 
required for the onset of the dynamo does not depend on the wall conductivity parameter 
c a . This is because c a drops out of the problem when 7 = 0, as follows from Eq. (20). 
However, c a affects the dynamo growth rate: larger values of c a (larger wall conductivity) 
lead to lower |7|. Therefore, in the limit of a perfectly conducting shell no growing field is 
possible, since 7 — >> 0. 

The wall permeability parameter c M has a strong influence on both dynamo threshold 
Rm cr and growth rate 7. A ferritic wall facilitates dynamo action. As shown in Fig. 4, the 
critical magnetic Reynolds number Rm cr decreases with increase of c M : from Rm cr ~ 244 
when c M = to Rm cr ^154 when c M — >> 00. These results are consistent with previous 
theoretical dynamo studies in other geometries [9, 10, 12, 13], which indicated reduction of 
Rm cr for the ferritic- wall boundary conditions. 

Estimates for typical parameters of MPDX show that its wall is very conducting and 
non- ferritic with c a « 30 and c M « 0. Under these conditions, dynamo action is achievable 
for the considered flow if Rm > 244, the respective dynamo growth rate at Rm = 400 is 
r « 3.6 s" 1 . 

In summary, we have studied the influence of finite conductivity and permeability of 
the wall on a plasma dynamo in a sphere. Our results show that the dynamo threshold is 
affected only by the wall permeability, while the dynamo growth rate depends on both wall 
properties. 

The authors wish to thank C. Sovinec for valuable help and discussions related to 
NIMROD. 
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